Manifestation of photonic band structure in small clusters of spherical particles 
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We study the formation of the photonic band structure in small clusters of dielectric spheres. 
The first signs of the band structure, an attribute of an infinite crystal, can appear for clusters of 5 
particles. Density of resonant states of a cluster of 32 spheres may exhibit a well defined structure 
similar to the density of electromagnetic states of the infinite photonic crystal. The resonant mode 
structure of finite-size aggregates is shown to be insensitive to random displacements of particles 
off the perfect lattice positions as large as half-radius of the particle. The results were obtained 
by an efficient numerical method, which relates the density of resonant states to the the scattering 
coefficients of the electromagnetic scattering problem. Generalized multisphere Mie (GMM) solution 
was used to obtain scattering matrix elements. These results are important to miniature photonic 
crystal design as well as understanding of light localization in dense random media. 
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The idea of employing photonic crystals, structures 
with periodically modulated dielectric constant, to ma- 
nipulate the density of states (DOS) of electromagnetic 
waves (EM) jl], |2| led to an explosion of interest, both 
academic and practical ||, §]. In a photonic crystal, in a 
given crystallographic direction, light can propagate only 
for the frequencies within photonic bands, described by 
a dispersion relation lo = w(k), similar to electron de 
Broglie waves in crystals. If the gaps between the bands 
overlap for all directions in the crystal, photonic crystal 
posses a complete photonic band gap || . For frequencies 
inside the complete gap, density of EM states turns to 
zero, that leads to new phenomena such as suppression of 
spontaneous emission, light localization, zero-threshold 
lasers, all-optical transistors and circuitry, and anoma- 
lous nonlinear properties jl], ^, ||, [|, |6|. For theoretical 
considerations, the size of the crystal is usually assumed 
infinite. In practice, however, one always deals with fi- 
nite structures. Naturally, a question arises: How large 
a photonic crystal should be, in order to exhibit photonic 
band structure? In this letter, we show that the modifi- 
cation of the density of EM states may begin for a clus- 
ter of 5 spheres. This suggests that scattering of light 
by ordered aggregates of spheres can be important in the 
problem of light propagation in dense random media. In- 
deed, an ensemble of mono-disperse, resonant scatterers 
(e.g. spheres) is considered the most favorable for 
satisfying Ioffe-Regel criterion for Anderson localization 
khcat ^5 1 1^7 @i H> H> 0' where k is the wavenumber 
in the medium and l sca t is the scattering length. Closely- 
packed face centered cubic (FCC) arrangement often ap- 
pear in self-assembled structures and it is shown to have 
lowest free energy. This fact points to a high probability 
of formation of ordered clusters in a collection of spheres. 
Even though the complete photonic gap can never be re- 
alized in an FCC structure of dielectric spheres in air, 
for some parameters, the total density of EM states can 
be significantly suppressed in certain frequency regions. 
Our results suggest that such suppression can occur due 
to the presence of the clusters even without long range 



order. 

Treatment of photonic structures, with all 3 dimen- 
sions of the order of the wavelength has been a challenge. 
When the distance between scatterers is comparable to 
their size, analytical methods, such as single scattering 
approximation, dipole approximation become invalid [Q. 
Numerical methods tailored for calculation of transmis- 
sion through a slab jL2| |l3) cannot be applied because 
it requires the lateral size of the slab to be larger than 
all characteristic lengths. Finite difference time domain 
(FDTD) method is used extensively in photonics design 
Q [l5j. Obtaining the density of states with FDTD- 
based ONYX method O] requires time-dependent solu- 
tion. Combined with dense A/20 spacial grids, it makes 
the problem computationally demanding in 3D, even for 
the smallest structure considered in this paper. However, 
the tendency to miniaturization of photonic devices calls 
for an efficient tool for dealing with such structures. An 
efficient alternative approach is needed. 

For small truncated photonic crystal structures, scat- 
tering language becomes appropriate. We use the GMM 
solution |ll| to obtain the elements of the scattering ma- 
trix of the collection of dielectric spheres. The density 
of resonant EM states of the cluster can be easily ex- 
pressed in terms of the scattering matrix. Within this 
approach, one only needs to evaluate vector-spherical- 
function expansion coefficients, there is no need, as in 
FDTD method, to find EM fields at every spacial point. 
The method requires considerably less memory and com- 
putation time resources, that makes it possible to use on 
a personal computer for up to 300 particles. An impor- 
tant class of photonic crystals, opals, is treated in this 
framework as an example. 

The main results of this paper are: (i) the signatures 
of the photonic band structure can appear for aggregates 
as small as 5 particles; (ii) in a cluster of 32 particles the 
density of (resonant) states can show strong resemblance 
to that of an infinite structure, with a pronounced de- 
pletion in the region of pseudo gap; (iii) photonic band 
structure of the finite clusters is tolerable to a certain 
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degree of random displacements of the particles off their 
lattice sites. 

We consider clusters of 5 and 32 spheres, with diame- 
ter 0.95 cm and refractive index n = 3.14. For both clus- 
ters we assumed FCC ordering with the distance between 
nearest neighbors of 1.9 cm. The system was chosen to 
match the setup used in the light localization experiment 
- Ref. m Smaller cluster is comprised of three planar lay- 
ers perpendicular to the 111 direction, with 1, 3, and 1 
sphere in each layer. In the larger cluster, particles were 
arranged in 5 layers with 3, 7, 12, 7, and 3 spheres in 
each layer, with the same stacking as in the first cluster. 

The cluster is an open system that does not sup- 
port stable modes, quasi-states of EM field are leaky, 
and density of the states cannot be specified. For elec- 
trons JltJ in open stochastic systems, Wigner delay time, 
t w {E) = h/i dln(S)/dE, expressed in terms of scattering 
matrix S, can be shown to be proportional to the density 
of resonant modes 
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where u>j and 7j are the energy and the linewidth of the 
jth resonance. For an infinite system, the poles of the 
scattering matrix S approach the real energy axis and the 
stationary modes are formed, jj become zero. For EM 
waves, one can define dwell time, r c ;, that closely follows 
EM version of the Wigner time, giving asymptotically 
the same result on the resonances . Dwell time can be 
readily found in terms of scattering coefficients as follows: 



VjV Ef 4 W{v^)dv 

c a scat {uj) 



(2) 



where W(r, u) is electromagnetic energy normalized to 1 
in vacuum, a sca t is the scattering cross-section, normal- 
ized by the total geometrical cross-section J7 7rrf, and 
c is speed of light. The integral in Eq. (||) is taken over 
the volume of each scatterer, Si. It gives in units of 
r/c and numerically coincides with cavity quality factor 
Q. The above expression can be understood physically as 
a ratio between the energy stored in a scatterer divided 
by outgoing energy current. J s W(r,u;)dr can be found 
by introducing an infinitesimally small absorption in the 
refractive index of the spheres n + 



/ W(r,uj) ~ —n lim 
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where Q a bs % is the absorption cross-section, and %i — 
1"KTij\ is size parameter of the ith sphere. The prob- 
lem of finding quasi-DOS of the finite cluster is now re- 
duced to finding the scattering coefficients Qabs z, Oscat- 
This step can be done by evaluating the expansion sc- 
ries of GMM solution ]l6| ]. To obtain a measure of the 
total density of resonant states r^ ot , expression Eq.(0) 
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FIG. 1: Photonic crystal band structure (middle) is com- 
pared to the dwell times obtained for scattering from a single 
sphere (bold solid line), cluster of 5 spheres (dashed line), and 
the cluster of 32 spheres (thin solid line). Upper graph cor- 
responds to EM plane wave incoming in 100 (X), the lower 
graph — 111 (L) crystallographic direction. 



was averaged over 9 symmetrically non-equivalent ori- 
entations (altogether 108 directions, accounting symme- 
try). We would like to stress, that within this frame- 
work, the source of error are numerical accuracy of the 
evaluation of expansion coefficients and limited number 
of orientations used in angular averaging, both of them 
can be controlled. 

One characteristic property of a photonic crystal is 
spacial dispersion, w(k). On Fig. |l| we compare dis- 
persion curves of the infinite photonic crystal, obtained 
using MIT photonic bands code fl8| , to the dwell times 
for two clusters, described above, calculated for 2 differ- 
ent incident angles corresponding to X and L crystallo- 
graphic directions. It can be easily seen, that even for 
smaller cluster of 5 particles shows the formation of 
the resonant modes at the positions of photonic modes in 
L-direction. For the cluster of 32 particles, well defined 
mode structure appears in both directions. However it is 
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more pronounced for L direction, where 6 resonant states 
between 12 and 16 GHz can be traced to L-modes of the 
infinite photonic crystal. The shift of the modes can as- 
cribed to the coupling between different modes owing to 
the finite size of the system. Indeed, in the infinite struc- 
ture the conservation of momentum of EM wave, due 
to periodicity, would decouple L-modes, whereas the fi- 
nite coupling in truncated crystal leads to the shift and 
widening of the resonances. Based on symmetry con- 
siderations, we can understand the stronger L-modes as 
compared to X-modes. The former are formed due to 
Bragg reflections from 3 or 5 (for the smaller and larger 
clusters respectively) planes formed by adjacent (nearest) 
neighbors in 111 direction, while in X direction, Bragg 
planes are formed only by next-nearest neighbors, which 
would make them more susceptible to the truncation. 

T^ ot and angularly (as for T^ ot ) averaged a scat in the 
region of the first two single sphere Mie resonances are 
shown on Figs. |2|(b,c). For comparison we included the 
total density of EM states, Fig. ^(a), calculated for the 
infinite crystal using the method described in Ref. |l^. 
Judging from scattering cross-section only, one cannot 
compare a finite cluster to the infinite structure, while 
T^ ot allows straightforward comparison with DOS of the 
photonic crystal 

Fig. ||(b) shows twofold decrease of a scat at Mie reso- 
nances in the cluster. This suppression of the scattering 
efficiency can be attributed to the hybridization of the 
single particle resonances ]2(J due to multiple scatter- 
ing. Dwell time is closely related to the near-fields [Q in- 
side the clusters, where formation of the Bragg standing 
modes leads to substantial modification of the spectrum 
- Fig. ||(c). For 32-sphere cluster one can see strong 
resemblance between the density of the resonant modes 
expressed by r^ ot and DOS in the photonic crystal, Fig. 
||(a). It is interesting to compare our results to ID case 
considered in Ref. ^l[ For a stack of periodically ar- 
ranged quantum wells it was argued that with increase of 
the system size, so-called, subradiant EM modes formed 
stable modes of the photonic passband with large life- 
times, similar to our result for dwell times in the 3D 
clusters. 

To assess the sensitivity of our results to disorder in 
the cluster, we performed the calculations of a sca t and 
T^ ot , Fig ||(d,e), for: the cluster of 32 spheres displaced 
by the half of the particle radius in random directions 
(dashed lines) off their initial positions and the cluster of 
32 spheres where 5 spheres replaced by the defect spheres 
of a different diameter (0.85 cm) - solid lines. One can 
see, that two types of disorder had different effect on the 
spectra. While positional (topological) disorder led only 
to smearing of only sharp resonances, the defects also in- 
troduced a number of new peaks. The later effect can 
be related to the new resonances introduced by the de- 
fect spheres. The effect of the topological disorder on 
photonic bandgaps was studied in (finite but comprised 
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FIG. 2: (a) Density of EM states in the infinite photonic crys- 
tal, (b) Normalized scattering cross-section for a single sphere 
(dotted), and clusters of 5 (dashed) and 32 (solid) spheres, (c) 
Normalized dwell time for the same as in (b). (d) Normalized 
scattering cross-section for cluster of 32 spheres, displaced 
by half-radius in random direction off their lattice position 
(dashed line), and a cluster of 32 spheres, where 5 spheres 
were replaced by the smaller, 0.85 cm diameter, spheres (solid 
line), (e) Delay time for the same as in (d). 
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of large number of particles) 2D £3| and 3D pho- 
tonic crystals Our conclusion on stability of the 
band structure to the topological disorder seems to be in 
line with Ref. |22|. However, besides obvious difference 
in the dimensionality of the space (we work in 3D), an- 
other important difference exists. In large systems, to 
explain the existence of the band structure one can re- 
sort to the periodic-on- average argument, while in such 
a small system as ours one cannot. 

It is worth noting, that FCC opal structures are not 
usually considered as a good candidate for potential ap- 
plications connected to complete photonic bandgap p5| . 
However, as it can be seen from our example, for some pa- 
rameters, DOS can be significantly suppressed in a wide 
spectral region. Fig. || also suggests that such structures 
can be tolerable to strong topological disorder. 

As it was already mentioned in the introduction, our 
results are relevant for the problem of light propagation 
in a dense random media made of particles with a narrow 
size dispersion. Topological disorder does not prevent 
the occurrence of ordered clusters. In Ref. ^, random 
displacement of scatterers off their lattice positions did 
not have a significant effect on the density of EM states. 
Our calculations performed for the clusters with similar 
random displacements (Fig. |2|(d,e)), confirm small effect 
on the formation of the photonic band structure. This 
also suggests that the light localization observed in Ref. 
^ may have been facilitated by depletion of the density 
of electromagnetic states due to the presence of the pseu- 
dogap. 

It can be further argued that when the scattering 
length l sca t becomes of the same order as the mean dis- 
tance between the particles, the occurrence of clusters 
cannot be ignored. Indeed, the scattering inside the clus- 
ters is dominated by Bragg reflection mechanism, while 
I Bragg can be as small as a few lattice constants for high 
index contrast structures. Furthermore, in Ref. [ p6| it 
was shown that periodicity of a photonic crystal may 
also exhibit itself in coherent back-scattering effect. The 
modification of coherent back-scattering due to the pres- 
ence of ordered clusters is an interesting problem that 
deserves a detailed consideration. 

In conclusion, we calculated scattering cross-section 
and density of resonant states in small 3D ordered aggre- 
gates of dielectric spheres within framework of the rigor- 
ous GMM solution. In contrast with previous studies we 
made no assumptions about the strength, size, or sepa- 
ration between scatterers. All multiple scattering effects 
were automatically retained in the solution. The results 
suggest that the photonic band structure of the infinite 
crystal can show up for clusters as small as 5 particles. 
Density of resonant states is significantly perturbed com- 
pared to single particle case. For a cluster of 32 particle 
the density of resonant states can show a close resem- 
blance to DOS of the infinite system. 
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